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A  method  is  devised  for  treating  the  transfer  of  a  light  particle 
such  as  a  hydrogen  atom  or  a  proton  between  two  heavy  particles, 
for  the  case  of  a  collinear  symmetric  reaction  AH  +  A  —  A  +  HA. 

Polar  coordinates  (p,  0)  with  origin  at  the  origin  of  the  usual 
mass- weighted  skewed- angle  potential  energy  contour  plot  are 
employed  to  simplify  the  calculation.  Symmetric  and  antisymmetric 
vibrationally  adiabatic  states  are  obtained  at  each  p.  They  provide 
an  effective  potential  for  the  p- motion  and  are  used  to  obtain  the  elastic 
adiabatic  phase  shifts  for  the  p-  motion  and  from  them  the  transfer 
probability.  The  p- motion  corresponds  approximately  to  that  of 
the  heavy  particles  and  the  0  to  that  of  the  light  one.  Semi- 
classical  limiting  expressions  are  obtained  and  provide  a 
description  of  light  particle  tunneling  at  low  energies  for  heavy- 
light-heavy  particle  systems.  Numerical  results  are  given 
and  the  present  approximation  is  good  for  the  systems  investigated. 
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FOREWORD 


The  reactive  transfer  of  a  light  particle  such  as  a  proton  or  hydrogen 
atom  is  of  interest  in  the  electrochemical  hydrogen  evolution  reaction  and 
in  other  systems.  In  the  present  report  an  approximate  method  is  for¬ 
mulated  for  treating  the  reaction  dynamics  of  transfer  of  a  light  particle. 

A  particular  intermolecular  potential  is  considered  to  permit  comparison 
with  exact  numerical  quantum  mechanical  results  for  the  transfer  probability. 

The  results  are  in  good  agreement  and  the  method  thereby  provides  a  new 
tunneling  description  for  the  light  particle.  The  present  report  is  restricted  to  the 

case  of  a  symmetric  reaction,  AH  +  A  —  A  +  HA.  The  work  is  being  extended 
to  asymmetric  reactions  AH  +  B  —  A  +  HB. 
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I.  INTRODUCTION 

In  low  energy  thermoneutral  or  nearly  thermoneutral  reactive  col¬ 
lisions,  only  one  vibrational  state  is  often  possible  in  both  reactants  and  products 
and  the  vibrational  quantum  number  appears  conserved.  For  these  (and  other) 
collisions  quasi- constants  of  the  motion  may  exist  which  remain  approximately 
conserved  throughout  the  collision.  If  such  a  quasi- constant  of  the  motion  can 
be  found  for  a  particular  reaction,  usually  in  the  form  of  a  suitably  defined 
generalized  vibrational  quantum  number,  it  can  improve  the  understanding  of 
the  dynamics  of  the  reaction  and  facilitate  its  computational  treatment. 

If  a  quasi- constant  of  the  motion  does  exist  for  a  particular  reaction  one 
of  the  ways  of  finding  it  is  to  find  a  set  of  coordinates  in  which  the  Hamiltonian 
is  nearly  locally  separable  and  one  of  the  coordinates  represents  a  generalized 
vibration  which  remains  to  a  large  extent  in  the  same  quantum  state  throughout 
the  collision. 

One  approach  of  this  kind  is  based  on  the  near- separability  of  the 
potential  energy  surfaces  for  some  reactions  in  the  vicinity  of  the  minimum 
energy  path  from  reactants  to  products.  In  terms  of  "natural  collision 
coordinates”,  *  designed  to  take  advantage  of  such  near  separability  of  the 
potential  energy  surface,  the  vibration  of  the  reactants  gradually  changes  its 
character  during  the  collision  to  become  a  vibration  of  the  products,  passing 
through  a  symmetric  stretch  vibration  or  the  collision  complex.  Such  a  near¬ 
separation  is  expected  to  be  most  suited,  other  things  being  equal,  to  col¬ 
lisions  in  which  the  minimum  potential  energy  path  has  a  low 
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curvature  in  the  usual  mass-weighted  Cartesian  coordinate  space.  It  has 

3 

been  useful  for  treating  a  variety  of  other  reactive  collisions. 
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One  class  of  reactions,  the  exchange  of  a  light  particle  between  two  heavy 

ones,  has  radically  different  dynamical  features  and  should  be  given  a  different 

treatment.  The  "reaction  path  curvature"  is  typically  large  in  some  regions 

because  of  the  acuteness  of  the  skewed  angle  in  the  usual  mass- weighted 

coordinate  plot  for  this  reaction  class.  The  large  difference  in  masses  in  this 

4 

case  points  to  an  effective  Born-Oppenheimer  type  adiabatic  separation  between 
the  motion  of  the  light  atom  relative  to  any  of  the  heavy  particles  and  the  relative 
motion  of  the  heavy  particles. 

The  adiabatic  separation  of  the  motion  of  heavy  and  light  atoms  has  been 

5 

successfully  utilized  in  finding  eigenvalues  of  hydrogen  bonded  systems.  In 
collision  problems  a  closely  related  approximation  of  treating  the  heavy 
particles'  motion  classically  and  the  light  particle  motion  quantum  mechanically 

C 

has  also  been  used.  A  semiclassical  treatment  which  utilizes  essentially  the 
same  dynamical  features  to  simplify  the  expressions  for  the  classical  S- matrix 

7 

has  also  recently  been  developed. 

In  the  present  paper  the  adiabatic  separability  of  the  light  atom  motion 
from  the  heavy  atom  motion  is  used  to  simplify  the  quantum  mechanical 
treatment  of  reactive  collisions  which  involve  exchange  of  a  light  particle, 
usually  a  hydrogen  atom  or  a  proton,  between  two  heavy  ones.  Only  the  sym¬ 
metric  exchange  of  hydrogen  atom  between  two  heavier  atoms  is  treated  in 
detail  in  the  present  paper.  Simple  approximate  expressions  for  the  reactive 
transition  probability  in  the  ground  vibrational  state  are  derived.  The  non- 
symmetric  case  will  be  discussed  in  a  later  paper. 

The  Schrfi dinger  equation  for  the  system  and  the  coordinate  system  in 
which  the  adiabatic  separation  can  be  utilized  are  presented  in 
Sec.  n.  The  transformation  of  the  SchrO dinger  equations  into  a  set  of  coupled 
ordinary  differential  equations  in  the  adiabatic  representation  is  given  in 
Sec.  m.  The  symmetric  exchange  case  and  the  adiabatic  approximation  which 
leads  to  a  simple  expression  for  the  reactive  transition  probabilities  are  given 
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in  Sec.  IV.  A  semiclassical  limit  expression,  valid  under  certain  conditions, 
is  derived  from  the  quantum  mechanical  result  in  Sec.  V.  Calculations  of  the 
transition  probabilities  for  two  model  potential  energy  surfaces  are  given  in 
Sec.  VI.  A  discussion  of  the  treatment,  its  possible  extensions,  and  a  classical 
path  approximation  are  given  in  Sec.  VII.  The  principal  equations  in 
the  paper  are  Eqs.  (4.10),  (5.7),  (5.11),  and  (7.  5). 


H.  SCHRODINGER  EQUATION  AND  THE  COORDINATE  SYSTEM 


Ths  Schrodinger  equation  for  a  collinear  reactive  atom-diatom  collision 
AB  +  C  —  A  +  EC  can  be  written,  after  separating  the  center  of  mass  motion, 
in  atomic  units  as: 


1 

2Ma 


2mc 


*  V(ra>  Ra) 


* 


Eip  ,  a  =  1,  2 


(2.1) 


where  r;  is  ZB  -  ZA,  (Z  denotes  position  along  line),  Ri  is  Zc  -  ZAB,  where 
Zab  is  the  position  of  the  center  of  mass  of  the  AB  molecule,  and  m,  and  M1 
are  the  corresponding  reduced  masses.  The  subscript  2  refers  to  the  same 
quantities  for  the  reverse  reaction.  So  r2  is  Zc  -  ZB  and  R2  is  ZBC  -  ZA, 

ZBC  being  the  position  of  the  center  of  mass  of  the  BC  molecule. 

For  systems  with  »  ma  (a  =  1,  2)  at  low  energy  a  Born-Oppenheimer^ 
type  adiabatic  separation  of  the  hydrogen  atom  motion  (along  ra)  and  the  heavy 
particle  motion  (along  R^)  may  exist  and  can  be  utilized  to  develop  an  approxi¬ 
mate  solution. 

In  order  to  simplify  the  derivation  it  is  convenient  to  mass-weight8 
the  coordinates  in  (2. 1)  as 

i  i 

xa  =  Ra;  ya  =  ma2  ra  (a  =  1,2)  (2.2) 
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The  resulting  typical  shape  of  a  potential  energy  function  for  a  collinear 
arrangement  of  two  heavy  atoms  with  a  light  one  between  them  plotted  in  mass- 
weighted  coordinates  is  shown  on  Fig.  1  as  a  contour  plot.  The  sharp  angle  y 
between  the  axes  on  the  plot  is  due  to  the  large  ratio  of  the  scaling  factors  in 
(2.  2),  i.  e. ,  to  a  large  mass  ratio  and  is  given  by 

1 

y  =  tan"1mB/(MQ,ma)2  ,  (2.3) 

where  m^  is  the  mass  of  the  middle  atom  B.  When  NT  /in  »  1  the  angle  y 

d  ot  at 

for  a  symmetric  reaction  is  approximately  (2mg/mA)2. 

The  Schro dinger  equation  in  coordinates  (2.  2)  becomes 


1 

1 


dx 


\  ~  +  V(x,y) 

2  3y 


4>  =  Eip  , 


(2.  4) 


It  is  the  same  as  the  SchrSdinger  equation  with  mass  1  moving 

in  a  two-dimensional- - 

potential  field  V(x,  y).  The  characteristic  feature  of  the  shape  of  the  potential 
energy  function  in  mass- weighted  coordinates  is  the  existence  of  two  long 
narrow  almost  parallel  channels  (Fig.  1).  The  motion  of  the  particle  along 
the  channel  length  corresponds  largely  to  the  relative  motion  of  the  two  heavy 
particles,  while  across  each  channel  it  corresponds  to  the  motion  of  the  light 
particle  relative  to  one  of  the  heavy  atoms. 

The  relevant  features  of  the  potential  energy  surface  in  coordinates 
(2.  2)  which  lead  to  vibrational  adiabaticity  can  be  illuminated  by  plotting 
an  example  of  a  low  energy  reactive  classical  trajectory  for  a  symmetric 
hydrogen  atom  exchange  reaction  shown  in  Fig.  1.  The  vibrational  state  of 
the  system,  i.  e. ,  the  classical  action  variable  for  the  motion  parallel  to  the 
channel  width,  remains  largely  unchanged  during  the  slow  approach  of 


k 
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the  reactants,  due  to  the  very  moderate  change  of  the  profile  of  the  bottom  of 
the  well  over  a  vibrational  period  of  the  trajectory.  Subsequently,  the  AB 
vibration  can  be  viewed  as  being  transformed  briefly  (for  a  half  a  vibrational 
period)  into  an  asymmetric  stretching  vibration  of  the  ABC  complex,  when  the 
trajectory  passes  from  the  reactants'  channel  into  the  products'  one.  Once 
in  the  product  valley  the  trajectory  remains  in  the  same  vibrational  state  during 
the  separation  of  the  heavy  atoms.  This  state  is  also  the  same  as  the  initial 
one  due  to  the  exact  matching  of  the  reactant  and  product  vibrational  potential 
wells. 


At  large  distances  between  the  colliding  partners  the  reactant  and  the 
product  vibrational  states  are  uncoupled  from  each  other  and  the  eigenfunctions  x 

of  the  vibrational  part  of  the  Hamiltonian  can  be  used  to  construct  the  asymptotic 
solution  to  the  Schrodinger  equation. 


1 

2 


d2 


dy 


+  V(xa,yflf) 


a 


xf  (ya)  -  <x“(ya)  (0-1. 2)  • 


(2.  5) 


The  value  of  a  =  1  refers  to  the  reactants'  channel  and  a  =  2  to  the  products'. 

At  sufficiently  large  Ri  and  R2,  one  has 

2 

*  =  T,  E  *?(xa)xf(ya)  •  (2-6) 

a=l  i  1  1  a 

where  <6?  (x  )  is  the  wavefunction  for  the  x  motion.  The  sum  over  i  is 
l  a  a 

over  all  asymptotic  states  in  channel  a. 

In  order  to  transform  the  Schrb dinger  equation  into  a  set  of  ordinary 
differential  equations  which  can  be  used  to  generate  solutions  of  type  (2.  6)  for 
large  xa  one  needs  a  pair  of  curvilinear  coordinates  which  can  be  made 
arbitrarily  close  to  xa,  ya  near  the  bottom  of  each  valley  of  the  potential 
energy  surface,  for  sufficiently  large  xa< 

Several  sets  of  coordinates  with  those  properties  have  been  proposed  in 

9-13  12  12 

the  literature.  In  this  paper  the  polar  coordinates  » 


f 
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P  =  (xa2+  y^)2 :  9  =  tan'1  (ya/*a)  (2.7) 

which  provide  a  combination  of  nearly  optimal  separability  and  extremely 
simple  form  of  the  Hamiltonian  will  be  used. 

For  large  p  and  moderate  ya  one  has 

P  »  xa  (2.8) 

and  the  vibrational  coordinate  y^  is  most  conveniently  expressed  in  terms  of 
the  arc  length  along  a  circular  arc  with  origin  at  x^  =  X^  =  0. 

s  ~  pd  ,  (p  —  oo)  f  (2.  9) 

namely  at  large  values  of  p  one  has 

y>  =  s;  y2  =  prs  .  (2. 10) 

(This  s  is  not  the  s  in  Ref.  1. ) 

III.  COUPLED  CHANNEL  EXPANSION 

Introducing  the  change  of  variable  (2.  7),  the  Schrodinger  equation  (2.  4) 
is  transformed  into: 


1 

1 


£_ 

dp2 


1  l  a 

2  p  3p 


+ 


H0(B,p) 


&(9,  P)  =  E  <h(B,p) 


(3.1) 


where  Ho(0,p)  is  the  quasivibrational  part  of  the  Hamiltonian  in  which  p  is 
a  parameter: 


H0 


1  _L  jL 

2  p2  d92 


+  V(9,p)  . 


(3.2) 


At  fixed  p  one  can  use  (2.  9)  and  transform  (3.  2)  to  the  form  of  a  one¬ 
dimensional  Hamiltonian  in  the  variable  s. 
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Ho  *  "  -2  -^T  +  V(s/p,p)  (3.3) 

3s 

H0Xj(s  ,P)  =  £j(p)Xj(s:p)  (3.4) 

where  the  x^sjp)  can  be  taken  as  real,  and  /  x^sjp)  Xj(s£)  p&G  equals  6. .  at 
fixed  p.  When  p  becomes  large,  we  shall  indicate  whether  s  is  in  the  range 
corresponding  to  the  reactants'  channel  by  writing 

Xj(s:p)  —  Xj^s-.p)  =  Xj‘(yi)  (P  — «)  ,  (3.5) 

where  x^i)  satisfies  (2.  5) ,  or  whether  it  is  in  the  range  corresponding  to  the 
products'  channel  by  writing 

Xj(s:p)  —  Xj2(s:p)  s  Xj2(y2)  •  (3-6) 

For  large  p  (3.2)  reduces  to  a  vibrational  Hamiltonian  for  a  diatomic 
molecule  in  the  vicinity  of  a  potential  minimum,  given  by  (2.5) ,  when  expressed 
in  terms  of  (x,  y)  coordinates. 

Since  the  solution  of  (3.2)  for  large  p  can  be  written  in  terms  of  the  solu¬ 
tions  of  (2.5) ,  one  can  seek  a  solution  of  (3.1)  in  a  form  which  can  be  readily 
reduced  to  the  form  of  the  asymptotic  solutions  of  (2.4)  for  large  p,  namely 

=  £  0,(p)  X^sp)  (3.  7) 

J  ]  J 

which  for  large  p  can  also  be  rewritten  instead  as  (2.  6) .  Wavefunctions  rj^t/jp) , 
defined  by 

iJ^Sp)  =  Ps  X^sp)  ,  (3.  8) 

with  the  normalization 

fylV  =  frfJjd*  =  ®ij  (3.9) 


are  next  introduced. 


s 


1 

Substituting  (3.7)  and  (3.8)  into  (3.1) ,  premultiplying  by  p1  ^(sp) ,  using 
(2.9)  and  integrating  over  6  one  gets  an  infinite  set  of  ordinary  coupled  dif¬ 
ferential  equations  in  p  for  the  coefficient  functions  <p^(p) 


1  d  ,  1  \ 


<f>^(p)  +  S  (-i  Py  ^  Qy)  <p  j(p)  =  0  ,  (3.  10) 


where  P--  =  <xt!-  *35  !xj> 


Qij  -  <X,|- 


dp 


2  ' 


X;)  • 


(3.  11) 
(3.  12) 


ei  is  given  by  Eq.  (3. 4)  and  the  brackets  denote  integration  over  6,  as  in  (3.  9) . 

Equation  (3.10)  is  formally  equivalent  to  the  Schrodinger  equation  (2.  4) 
or  (3. 1).  By  solving  it  and  imposing  the  scattering  boundary  conditions 
one  can  get  the  full  scattering  solution  of  (2.  4).  As  shown  in  the  next  section 
however,  for  symmetric  hydrogen  atom  exchange  between  two  heavy  atoms 
one  can,  at  low  energies,  obtain  a  good  approximate  solution  without  solving 
coupled  equations. 


IV.  SYMMETRIC  HYDROGEN  ATOM  EXCHANGE 

If  atoms  A  and  C  are  the  same,  the  potential  energy  function  V(p,  9)  is 
symmetric  about  6  =  y/2,  and  there  are  two  distinct  kinds  of  eigenfunctions 
of  (3.  3)  which  are  symmetric  and  antisymmetric  about  9  =  y/2,  respectively. 
All  the  coupling  elements  (3.  9)  between  the  symmetric  and  the  anti¬ 
symmetric  eigenfunctions  x^  of  (3.  3)  vanish  by  symmetry  and  the  system  (3.  8) 
can  be  decoupled  into  two  systems  for  the  coefficients  d>s(p)  and  <£a(p)  of  the 
symmetric  and  the  antisymmetric  wavefunctions  Xj(s,p),  respectively. 
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1  ^-,E-£f+-L) 

2  dp  1  8p 


W  -  jpQ  (pij  ®  *  2  «|>  *,«  =  * 


@  =  s,  a) 


(4.1) 


Further,  for  collisions  involving  transfer  of  a  hydrogen  atom  between  two 

heavy  atoms  at  low  energies  the  Hamiltonian  is  nearly  adiabatically  separable 

0 

in  the  coordinates  p,  9  and  the  coupling  elements  P0j  between  the  ground  and 

the  higher  states  of  the  same  symmetry  are  expected  to  be  small.  In  addition, 

at  low  energies  near  the  reaction  threshold  only  one  state  of  each  symmetry  is 

open  asymptotically,  and  the  effective  coupling  P^.  to  the  classically  forbidden 

14 

states  is  weak.  Since  the  depend  solely  on  the  P^  the  infinite  sum  in 

(4. 1)  can  be  neglected  and  the  equation  reduces  to 

<6?  =  E<f>?  (a  =  s,  a)  ,  (4.2) 

where  <&  refers  to  p- motion  for  the  lowest  vibrational  state  of  symmetry  a. 
The  solution  of  (4.  2)  for  large  p  yields  the  scattering  solution  of 

(3. 1)  in  terms  of  the  symmetric  and  antisymmetric  ground  state  eigen¬ 
functions  of  (3.  3).  The  latter  can  be  written  as  linear  combination  of  the 
ground  vibrational  states  of  the  reactants  and  products. 

xf,a(s:p)  =  —  [Xo(s;P)f  Xo  (s:p)  ]  .  (4.3) 

n 

The  Xo  (s.P)  corresponds  to  the  system  being  in  the  reactants',  and  i&isw)  in 
the  products'  configuration.  The  s  superscript  in  (4.  3)  refers  to  the  plus  sign. 

For  sufficiently  large  p  the  coordinates  (s,  p)  can  be  made  arbitrarily 
close  to  the  Cartesian  coordinates  xa,  ya  in  the  region  near  the  minimum  of 


1 

1 


dp 


+  € 


a 


8p 


11 


the  well  in  the  reactant  (a  =  1)  or  product  (a  =  2)  configuration.  The  vibra- 

1  2 

tional  wavefunctions  x^  ’  (s;p)  then  reduce  to  the  uncoupled  vibrational  wave- 

functions  of  the  reactants  and  the  products,  as  in  the  second  half  of  (3.  5) 
and  (3.  6). 

The  wavefunction  solution  of  (3. 1)  in  the  ground  vibrational  state 
!^(s,p)  =  <f>f(p)  xf(s:p)  +  (bo  (p)  Xo(s:p)  (4.4) 

can  hence  be  written  for  large  p,  using  (4.  3),  as 

>p(s,p)  ~  (t>o(p)  Xo(s-P)  +  <t>20(p)  Xo(s:p) 

~  <2>i(xt)  Xo(yi)  +  </>o(x 2)Xo(y2)  •  (4.5) 


Comparing  (4.  3),  (4.  5),  (3.  5)  and  (3.  6)  one  can  write  for  large  p 

0o’2(Xi>2)  ~  |$J(p)  ±  d>o  (P)J  , 


(4.6) 


where  the  1  superscript  refers  to  the  plus  sign. 

Writing  the  asymptotic  solution  of  (3. 1)  in  terms  of  an  S-matrix 


*  =  (e-*xi 


(4.7) 


where  S^0  and  are  the  amplitudes  of  the  nonreactive  and  reactive  outgoing 
waves,  respectively,  i.  e. ,  the  nonreactive  and  reactive  S-matrix  elements,  one 
obtains 


=  l  (e2i*°  +  e2i*°)  =  ei(*«  +  *»  *  cos(£f  +  £) 


=  I  <e2i^°  -  e2i*°)  =  ei^®  +  5®)sin(«f  -  Iq)  , 


eta 


s  ,  ea, 


(4.  8) 
(4.  9) 


where  ,  lo  are  the  elastic  phase  shifts  obtained  by  solving  (4.  2). 
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The  zero  order  approximation  to  the  reactive  transition  probability  in 
a  symmetric  collinear  hydrogen  atom  exchange  collision  is  hence  given  by 

PS  =  |sSf2=  sin2(|f  -  *a)  .  (4.10) 

Equation  (4. 10)  is  the  same  as  the  expression  for  the  S-matrix  elements  in 
case  of  exact  resonance  in  atomic  collisions  involving  transfer  of  an  electron 
or  electronic  excitation  between  two  identical  atoms.  ^ 

Equation  (4. 10)  is  the  desired  expression  for  the  probability  of  trans¬ 
ferring  a  light  particle  in  a  low  energy  collinear  collision  based  on  the 
adiabaticity  of  the  light  particle  vibration  throughout  the  collision.  In  order 
to  obtain  numerically  the  values  of  the  transition  probability  as  given  by 
(4. 10)  as  a  function  of  energy  at  low  energies  one  needs  to  solve  first  for  the 
two  lowest  eigenvalues  ef’  a(p)  of  H0,  Eq.  (3.  3),  for  several  values  of  the 
parameters  p  and  then  evaluate  the  phase  shifts  ’ a  for  scattering  on  the 
effective  potentials  ef ’ a  for  the  heavy  particle  motion  along  p  by  solving 
Eq.  (4.  2)  at  the  desired  energies. 

Some  numerical  results  of  applying  Eq.  (4. 10)  are  given  in  Sec.  VI. 

In  the  next  section  we  evaluate  the  semiclassical  limit  of  the  quantum 
mechanical  result  (4. 10),  thereby  exhibiting  some  of  the  relevant  aspects  of 
the  collision  dynamics  in  a  more  intuitive,  semiclassical  way.  (Cf  also  Sec.  VII. ) 

V.  SEMICLASSICAL  LIMIT  OF  THE  TRANSITION  PROBABILITIES  EXPRESSION 

If  the  potential  energy  surface  for  the  collinear  symmetric  hydrogen 
exchange  V(s:p)  has  a  sufficiently  large  barrier  the  reaction  will  at  low 
energies  proceed  by  tunneling  and  the  probability  for  reaction  at  energies 
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below  the  classical  threshold  will  be  nonzero  but  small.  All  the  quantities 
needed  for  calculating  the  transition  probability  in  (4.  10)  for  those  energies 
can  be  evaluated  semiclassically. 

The  "vibrational"  potential  energy  curve  V(s;p)  will  be  a  double  minimum 
potential  as  a  function  of  s  for  all  the  values  of  p  relevant  for  determining  the 
reaction  probabilities  in  the  tunneling  region.  Moreover,  both  the  lowest  sym¬ 
metric  and  the  lowest  antisymmetric  eigenvalue  of  the  double- well  potential 
will  be  below  the  barrier  maximum  (Fig.  2)  so  that  a  simple  semiclassical 
approximation  can  be  used  for  their  evaluation.  The  semiclassical  condition 

for  determining  the  lowest  pair  of  eigenvalues  of  a  symmetric  double-well 

1  fi  17 

potential  below  the  top  of  the  barrier  is  given  by1  ’ 

f  2  p(s;p)  ds  =  »/ 2  t  (5.1) 

si 

where 

s3  i 

Q(p)  =  /  p(s;p)  ds,  p(s;p)  =  j  2  [e  (p )  -  V(s;p)]|3  .  (5.2) 

sz 

The  Sj  denote  the  turning  points  for  the  motion  on  the  double- well  potential 
(Fig.  2)  and  are  slightly  different  for  the  +  and  -  states  in  (5.  1) .  The  s^ 
denote  the  turning  points  for  the  motion  in  each  of  the  two  separate  single- well 
potentials  (Fig.  2) .  The  mean  energy  of  the  first  two  eigenvalues,  which  will  be 
denoted  by  e0(p)  is  in  the  limit  of  small  splitting  given  by  the  usual  semiclassical 

eigenvalue  condition  for  a  single  well  potential1® 
o 

f  2  p(s^p)  ds  =  jr/2  (5.3) 

-s° 

16  17 

Their  splitting  is  given  by  ’ 

4l-f!W=fE.W«',W  •  (5-4) 


Here,  E0(p)  is  the  ground  state  eigenvalue  for  a  single  well  potential  as 
determined  by  (5.  3) .  E0(p)  is  measured  from  the  bottom  of  the  well,  i.  e. , 


E0  (p)  =  €0  (p)  -  V'min(s:p).  The  energy  splitting  of  the  eigenvalues  of  sym¬ 
metric  double  well  potential  is  hence  proportional  to  the  tunneling  amplitude 
from  one  well  to  the  other  exp[-Q(p)  ]. 

If  the  barrier  on  the  potential  energy  surface  is  sufficiently  large 
both  the  symmetric  and  the  antisymmetric  eigenvalues  e  ’  (p)  will  change 
monotonically  with  increasing  p  and  the  elastic  phase  shifts  for  scattering  on 
€  (p)  and  e  (p)  can  be  evaluated  semiclassically  as 


^  =  /<p)p  -  f  pV)dp':  0  =  s, a)  ,  (5.5) 

p3 

where  p^(p)  =  {  2[  E  -  e^(p)]pand  p^  is  the  classical  turning  point  on  the 
effective  potential  c^(p)  for  the  p  motion.  The  difference  of  the  elastic 
phase  shifts,  which  determines  the  reactive  transition  probability  in  (4. 10), 
is  then 


£S  -  £a  =  /  ps(p)dp  -  f  pa(p)dp;  p  —  ® 


(5.6) 


A  difference  of  phase  integrals  similar  to  (5.  6)  appears  in  the  semi- 

1  9 

classical  treatment  of  curve  crossing  problems  and  is  often  approximated  by, 


(o-  (f  =  f  (5.7) 

i 

where  p0(p)  -  {2[E-  €0(p)]}  5  and  p0  is  the  value  of  p  at  the  classical  turning 

point  of  the  ’’mean"  potential  energy  curve  €0  (p).  Equation  (5. 7)  is  often 
19 

referred  to  as  the  "classical  path"  approximation  since  it  leads  to  results 
similar  to  those  obtained  by  assuming  the  existence  of  the  time  dependent 
classical  path  p(t).  The  assumptions  used  in  obtaining  (5.7)  are  however 


fi. 


i 
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less  stringent  than  the  assumption  of  the  existence  of  a  time-dependent 

19c 

classical  path  for  the  p  motion. 

Using  (5.  4)  and  (5. 10)  can  be  written  as 

r>R  _  lt$  _  /fs  _  ta\2 


where 


PS  =  sin2(4o  -  O  *  (4o  -  O 


s  a  .2  r  W  -Q(p)  . 

'  4°  '  *  ■>„.  S5T  *  • 


(5.  8) 


(5.  9) 


Equation  (5.  9)  gives  the  semiclassical  limit  of  the  hydrogen  atom 
exchange  probability  as  an  integral  of  the  hydrogen  atom  tunneling  amplitude 
over  the  configuration  space  of  the  heavy  atoms.  A  further  physical  insight 
into  the  structure  of  (5.  9)  can  be  obtained  by  the  explicit  introduction  of  the 
classical  path  approximation  as  shown  in  Sec.  VTL 

We  next  proceed  to  evaluate  the  integral  in  (5.  9)  approximately.  Typically 
e“Q(p)  will  decay  with  increasing  p  sufficiently  rapidly  so  that  the  integral 
above  is  determined  by  the  values  of  the  integrand  in  a  narrow  region  around 
the  classical  turning  point  for  the  heavy  particles  motion  p0.  In  the  vicinity  of 
p0  the  slowly  varying  vibrational  eigenvalue  E0(p)  can  be  approximated  by  its 
value  at  the  turning  point,  E°.  Both  p*(p)  =  2[E-e0(p)]  and  Q(p)  are  next 
approximated  by  their  power  series  expansions  around  p,  truncated  at  the 
linear  term.  Thus, 


EoW  *  3 


(5. 10a) 


P o  (p)  *  2(p-p0)€0'  ;  €q  =  - 


de0(p) 

dp 


(5. 10b) 


Q <P)  *  Qo  +  (P-Po)Qo  :  Qo  = 


(5.  10c) 


P=Po 


where  Q,>  denotes  Q(pq)  .  Substituting  (5. 10)  into  (5.  9)  and  performing  the 
integration  one  obtains 
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rtE°0 


e~Q>  =  ye'Qo 


(5.  11) 


where  y  is  defined  by  the  right-hand  equality  in  (5.  11),  Qq  is  given  by 
(5.  2)  evaluated  atp  =  p0  and  ^(p)  is  defined  as  the  eigenvalue  satisfying  (5.  3) . 
Equation  (5. 11)  is  given  also  in  Ref.  7,  with  a  slightly  different  tunneling  path 
and  €q(p)  function  (polar  coordinates  were  not  used) . 

The  reactive  transition  amplitude  is  hence  given  by  the  penetration 
amplitude  for  the  hydrogen  atom  tunneling  e"^°  at  the  classical  turning  point  of 
the  heavy  particles  multiplied  by  a  preexponential  factor  y. 


VI.  NUMERICAL  TESTS 

The  reactive  transition  probabilities  as  a  function  of  the  collision 

energy  were  evaluated  for  two  different  LEPS  potential  energy  surfaces.  The 

20a 

first  one  was  chosen  using  as  asymptotic  parameters  those  for  the  F  +  HF  — 

FH  +  F  reaction.  Three  different  Sato  parameters  0.15,  0.1  and  0  were  used 

to  obtain  surfaces  with  a  low  (~  1.5  kcal) ,  moderate  (~6  kcal)  and  high  (~19  kcal) 

barrier  at  the  saddle  point.  These  three  surfaces  will  be  denoted  as  surfaces 

la,  lb  and  lc,  respectively.  Only  the  last  one  has  a  barrier  height  which 

roughly  corresponds  to  the  expected  barrier  for  the  reaction;^®*5  the  other  two 

were  included  to  examine  the  influence  of  the  barrier  height  on  the  reaction 

dynamics.  The  second  surface  has  parameters  chosen  to  simulate  the 

CH«  +  CH3  —  CHj  +  CH4  reaction.  The  asymptotic  parameters  are  the  same  as 

21 

those  of  Ref.  7a  and  the  Sato  parameter  was  chosen  to  be  s  =  0.  185.  The 

resulting  barrier  at  the  saddle- point  is  13.  8  kcal. 

The  quasi- vibrational  eigenvalues,  L  e. ,  the  eigenvalues  of  (3.  3)  were 

22 

evaluated  using  the  finite  element  method  of  Malik,  Eccles  and  Secrest 
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The  first  two  symmetric  and  antisymmetric  state  eigenvalues  are  shown  in 
Fig.  3  as  a  function  of  p  for  surface  lb. 

The  present  approximation  to  the  reactive  transition  probability 

23 

(4.  10)  was  obtained  by  solving  (4.  2)  numerically  using  Gordon’s  method 
to  obtain  the  elastic  phase  shifts.  The  reactive  transition  probability  as 
a  function  of  the  total  energy  is  shown  in  Fig.  4  for  surface  la,  b,  c  and  in 
Fig.  5  for  surface  2.  This  transition  probability  as  a  function  of  energy 
is  shown  on  a  logarithmic  scale  in  the  tunneling  region  in  Figs.  6  and  7. 
Some  classical  trajectory  results  for  the  moderate  barrier  surface  lb  are 
also  shown  in  Fig.  4. 

A  method  well-suited  for  performing  accurate  quantum  mechanical 
calculations  involving  a  symmetric  hydrogen  atom  exchange  between  two 
heavy  atoms  (or  any  atom  diatomic  molecule  reaction)  has  been  developed 
by  Kuppermann,  Kaye  and  Dwyer.  ’  The  results  of  the  present  approxi¬ 
mation  for  the  systems  studied  here  are  in  good  agreement  with  preliminary 

13c 

results  obtained  using  that  method.  In  particular,  the  shape  of  the 

transition  probability  vs.  energy  curve  for  all  the  systems  studied  is 

reproduced  quite  accurately.  The  features  included  in  this  curve  are  the 

26  27 

peaks  for  surface  la,  b  due  to  shape  resonances  in  the  solution  of  (4.  2)  ’ 

and  the  large  oscillations  (Figs.  4c,  5)  due  to  smooth  increases  of  the  phase 
difference  in  (4.  10)  through  2v.  There  is  a  small  shift  of  0.  04  -  0.  12 
kcal/mole  of  the  approximate  curve  in  Figs.  5-8  toward  the  lower 
energy  relative  to  the  coupled  channel - 
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results  in  the  threshold  region.  The  shift  increases  at  higher  energies, 
particularly  for  the  low  barrier  reactions.  Even  smaller  shifts  can  be 
expected  for  collisions  in  which  a  hydrogen  atom  is  transferred  between  still 
heavier  atoms.  In  contrast,  for  the  collisions  of  the  X  +  HX  series,  where 
X  is  a  hydrogen  atom  isotope  for  which  accurate  coupled  channel  results 

i  n  _  j 

exist, iJa’  a  the  corresponding  shifts  are  much  larger  (0.  5  to  1.  0  kcal  in  the 
24 

threshold  region) . 

VTI.  DISCUSSION 

The  results  presented  in  the  preceeding  section  indicate  that  there  is 
an  effective  adiabatic  separation  between  the  fast  hydrogen  atom  vibration  and 
the  slow  translation  of  heavier  atoms  in  hydrogen  atom  exchange  reactions  at 
energies  near  the  reaction  threshold.  The  polar  coordinates  used  here  are 
sufficiently  close  to  the  dynamical  variables  between  which  the  optimum 
adiabatic  separation  exists  at  any  point  in  the  reaction,  and  so  are  useful 
for  these  approximate  calculations. 

At  energies  at  which  only  the  ground  vibrational  state  is  asymptotically 

open,  the  reaction  transition  probability  for  a  symmetric  exchange  reaction 

can  always  be  written  rigorously  in  the  form  (4.  10)  if  one  replaces  the 

adiabatic  phase  shifts  a  obtained  by  solving  Eq.  (4.  3)  by  the  exact  phase 
25 

shifts  which  in  the  present  formulation  can  be  obtained  by  solving  (4.  2) . 
However,  only  for  the  transfer  of  a  light  particle  between  two  heavy  ones  are 
the  polar  coordinates  sufficiently  close  to  the  optimally  separable  ones  and 
(4. 10)  then  provides  a  good  approximation  to  the  transition  probabilities. 

The  principal  approximations  in  the  analysis  leading  to  Eq.  (4. 10) 
are  (i)  neglect  of  nonadiabatic  effects  and  (ii)  use  of  polar  coordinates  to 
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treat  the  dynamics.  Each  could  be  improved,  e.  g. ,  in  (i)  by  including 
perturbatively  a  nonadiabatic  correction  to  an  elastic  phase  shift  in  (4.  2), 
particularly  in  the  case  of  any  near  avoided  crossing  of  the  e0  (p)  and  e,  (p) 
curves.  However,  for  the  present  mass  ratios  the  current  approximations 
are  seen  to  be  quite  adequate  for  the  systems  treated. 

A  distinct  feature  of  the  hydrogen  atom  or  proton  transfer  reaction 
between  two  heavy  atoms  on  potential  energy  surfaces  with  high  barriers 
is,  as  shown  in  Figs.  5c  and  6,  the  unusually  large  amount  of  tunneling,  i.  e. , 
there  is  an  appreciable  reaction  probability  when  the  total  energy  of  the  reaction 
is  below  the  top  of  the  saddle  point  barrier.  For  reactions  with  different  mass 
combinations  for  which  the  characteristic  angle  y  in  Eq.  (2.  3)  is  much  closer 
to  90°  can  be,  at  least  qualitatively,  understood  in  terms  of  natural  collision 
coordinates.1  The  vibrational  motion  in  those  coordinates  is  perpendicular  to  the 
reaction  coordinate  and  the  part  of  energy  tied  up  as  vibration  is  unavailable 
for  crossing  the  barrier.  In  the  present  case  the  vibrational  coordinate  is 
parallel  to  the  reaction  coordinate  (see  Fig.  1)  only  in  the  region  of 
closest  approach,  and  is  utilized  in  surmounting  the  barrier  between  the 
reactant  and  the  product  channel. 

2ft 

We  now  turn  to  the  semiclassical  expression  (5.  9) .  This 
expression  points  to  an  intuitive  way  of  understanding  the  physical  factors 
which  determine  the  transition  probability  in  the  tunneling  region. 

The  value  of  the  transition  probability  in  the  tunneling  region  in  (5.  11) 
is  determined  mainly  by  the  amplitude  for  the  hydrogen  atom  tunneling  in 
one  dimension  across  the  ridge  between  the  two  channels  at  the  classical 
turning  point  for  the  heavy  particle  motion,  given  by  e"^°.  The  preexponential 
factor  y  is,  for  the  systems  lb,  lc  and  2,  within  a  factor  of  five 
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from  unity  and  depends  on  the  characteristics  of  a  particular  surface. 

Equation  (5. 11)  can  also  be  obtained  by  treating  the  heavy  particle 
motions  classically  and  the  hydrogen  atom  motion  quantum  mechanically. 

The  approach  used  above,  however,  yield  the  same  result  without 
the  restrictive  assumption  of  existence  of  a  time- dependent  classical 
trajectory  for  the  heavy  particle  motion.  The  classical 
path  picture  can  nevertheless  be  used  to  gain  additional 
insight  into  the  structure  of  (5.9)  and  into  the  nature  of  the  preexponential 
factor  in  (5. 11).  If,  starting  from  (5. 9), one  introduces  explicitly  a  classical 
trajectory  p(t)  for  the  classical  motion  of  the  heavy  atoms  on  the  one¬ 
dimensional  mean  potential  curve  e0  (p)  the  local  heavy  particle  velocity  and 
momentum  can  be  written  as 

Po  Co)  =  V0(p)  =  dp/dt  (m  =  1)  .  (7.1) 

The  local  vibrational  energy  E0  (p),  can  be  written  in  terms  of  the  vibrational 
frequency  j  hi/,  i.  e. ,  in  terms  of  the  vibrational  period  as 

E0fp(t)l  =  7r/r[p(t)]  (Ti  =  1)  ,  (7.2) 

where  r  is  the  vibrational  period  of  the  s  motion  for  a  given  p(t)  and  hence 
for  a  given  t.  Substituting  (7. 1)  and  (7.  2)  into  (5.  9)  one  obtains 

Is2l  =  f  e"Q(t)  d t/r  .  (7.3) 

-GO 

(A  factor  of  2  in  (5.  9) ,  and  symmetry  about  the  time  of  closest  approach 
[Q(t)  *  Q(-t)  ]  results  in  the  integration  limits  being  -*>  to  +«•  instead  of 
to  0. ) 

According  to  (7.  3)  the  reactive  transition  amplitude  can  be  obtained 
by  integrating  over  the  time-dependent  tunneling  amplitude  for  tunneling  of 
the  hydrogen  atom  into  the  product  region  over  the  time  of  approach  and 
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separation  of  the  heavy  particles. 

Equation  (7.  3)  is  a  classical  path  approximation.  However,  it  can  be 
used  to  obtain  a  semiclassical  tunneling  correction  for  use  of  an  ensemble  of 
purely  classical  trajectories  below  the  classical  threshold  as  follows. 
Dividing  the  integration  range  in  (7.  3)  into  increments  over  one  vibrational 


period  each  one  obtains 

I*  -  s  fYTie-M 
i  h 


(7.  4) 


Integrating  d (t/r)  over  an  interval  t.  +  r-  in  Eq.  (7.  4)  corresponds  in  an 
ensemble  of  trajectories  to  averaging  over  the  initial  vibrational  phase 
(classical  angle  variable) .  Thereby  Eq.  (7.  4)  becomes 


(7.  5) 


where  the  brackets  ( )  indicate  an  average  over  the  initial  vibrational  phase 
and  denotes  the  value  of  Q(t)  at  the  i'th  inner  vibrational  turning  point 
(i.  e. ,  the  point  closest  to  the  products'  region) ,  whose  p-value  in  turn  varies 
with  initial  vibrational  phase. 

Hence,  the  reactive  transition  amplitude  can  be  obtained  from  an 
ensemble  of  purely  classical  trajectories  by  summing  the  amplitudes  for 
tunneling  into  the  product  region  at  every  inner  vibrational  turning  point  during 
the  approach  of  the  reactants  and  during  their  separation. 

In  spite  of  the  additional  approximations  used  in  deriving  Eq.  (7.  5) , 
the  latter  could  give  better  results  for  the  transition  probability  in  the 
tunneling  region  than  the  quantum  expression  (4.  10.  The  use  of  exact  classical 
mechanics  provides  a  more  accurate  description  of  thep-s  dynamics  within  a 
channel  than  the  decoupling  approximation,  clearly  at  the  price  of  the 
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additional  effort  needed  to  solve  for  the  classical  dynamics  exactly. 

Another  advantage  of  Eq.  {7.  5)  is  that  its  use  can  straightforwardly 
albeit  heuristically  be  extended  to  multidimensional  systems  in  order  to 
obtain  an  approximation  for  the  tunneling  corrections  to  the  results  of 
purely  classical  trajectory  studies.  The  increasing  complications  with 
systems  with  many  degrees  of  freedom  make  the  exact  quantum  calculations 
computationally  difficult  at  the  present  time,  while  the  classical  trajectory 
studies  can  be  done  more  easily. 
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FIGURE  CAPTIONS 


Figure  1.  A  potential  energy  contour  plot  of  the  surface  lb  in  mass- weighted 
coordinates.  Several  values  of  the  contours  are  listed  in  kcal.  A 
reactive  classical  trajectory  with  energy  8.  93  kcal,  slightly  above 
the  classical  threshold,  is  also  shown  (dashed  line) .  The  sudden 
jump  of  the  trajectory  from  the  reactants'  to  the  products’  region 
well  before  reaching  the  saddle  point  is  characteristic  for  transfer 
of  a  light  particle  between  two  heavy  ones.  This  is  in  contrast  to 
the  usual  picture  of  a  vibrationally  adiabatic  collision  in  which  the 
trajectory  winds  along  the  minimum  potential  energy  path.  Most 
trajectories  are  less  symmetric  than  the  one  shown. 

Figure  2.  A  cut  through  the  potential  energy  function  V(s,p)  forp  =  13.7  on 
surface  2;  semiclassical  (WKB)  single-well  eigenvalue  ^(p) ,  as 
determined  by  Eq.  (5.  3)  (full  line) ;  semiclassical  approximation  to 
the  symmetric  and  the  antisymmetric  double  minimum  eigenvalues, 
given  by  Eq.  (5.  1)  (dashed  lines) .  The  points  si  and  s?  are  defined 
in  the  text. 

Figure  3.  The  first  two  symmetric  and  antisymmetric  diagonal  matrix  elements 
e^(p)  (Eq.  4.  3)  for  surface  lb  (moderate  barrier  surface  with 
F  +  HF  asymptotic  parameters) . 

p 

Figure  4.  The  ground  vibrational  state  reactive  transition  probability 

(Eq.  4.10)  (full  line)  as  a  function  of  total  energy  for  surfaces  la,  b,  c 
(FH  +  F  with  low,  moderate  and  high  saddle-point  barrier) .  The 
vertical  arrows  on  the  abscissa  denote  the  saddle-point  barrier 
heights  for  surfaces  b  and  c.  The  a  curve  has  two  sharp  peaks  as 
indicated.  Some  purely  classical  trajectory  results  for  the  moderate 
barrier  surface  lb  in  the  threshold  region  are  also  shown  (dashed  line) . 


Figure  5.  The  ground  vibrational  state  reactive  transition  probability  P00 
for  surface  2  (CH3  +  CH4  —  CH«  +  CH3) .  The  vertical  arrow  on 
the  abscissa  denotes  the  saddle- point  barrier  height. 

Figure  6.  Logarithm  of  the  reactive  transition  probability  P^,  vs.  total 
energy  (Eq.  4. 10)  for  surface  lb. 

p 

Figure  7.  Logarithm  of  the  transition  probability  P^,  vs.  total  energy 
(Eq.  4. 10)  for  surface  2. 
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